Commensurability effects for fermionic atoms trapped in ID optical lattices 
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Fermionic atoms in two different hyperfine states confined in optical lattices show strong com- 
mensurability effects due to the interplay between the atomic density wave (ADW) ordering and the 
lattice potential. We show that spatially separated regions of commensurable and incommensurable 
phases can coexist. The commensurability between the harmonic trap and the lattice sites can be 
used to control the amplitude of the atomic density waves in the central region of the trap. 
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Ultracold Bose or Fermi gases can be confined in ar- 
tificial optical lattices created by standing-wave laser 
fields. The low-energy properties of these systems can be 
described using models borrowed from condensed mat- 
ter systems [l| whose parameters and dimensionality can 
be controlled with high precision. This tunability has 
opened new avenues for understanding the physics of 
strongly correlated systems. Greiner et al. observed a 
superfluid to Mott insulator transition in a 3D optical lat- 
tice with bosonic ^"^ Rh atoms. 2] Interesting experimental 
results have also been obtained for fermions.[3:j 

The physics of cold fermionic atoms in optical lat- 
tices is predicted to be the one of the attractive Hub- 
bard model with the different hyperfine states playing 
the role of the spin states. In the case of one dimen- 
sion, bosonization predicts the formation of a Luther- 
Emery liquid for attractive interactions. Q The spin sec- 
tor is gapped, inducing exponential decay of spin correla- 
tions in contrast to singlet superconducting and charge- 
density wave correlations that have a power law decay. [5j 
Atoms are trapped in experiments and the interplay be- 
tween correlations, lattice, and confinement has to be 
properly addressed. Moreover, the study of the inho- 
mogeneity can be very relevant for the understanding of 
unconventional superconductors. For example, the den- 
sity of states of high-Tc superconductors presents spatial 
inhomogeneities.p Confined fermions in optical lattices 
are ideal experimental candidates to study these issues. 

Recently, Gao Xianlong et al. studied the unpolar- 
ized attractive Hubbard model in a ID optical lattice 
with harmonic confinement. Using the density matrix 
renormalization group (DMRG) approach, they observed 
the coexistence of spin pairing with an ADW that could 
be interpreted as a signature of a Luther-Emery liquid 
phase. This interpretation is further justified by the 
fact that the momentum of the ADW is proportional 
to the average atomic density in the bulk of the trap 
kADW — T^n. These ADWs could be detected measuring 
the elastic light-scattering diffraction pattern (the Fraun- 
hofer structure factor), proportional to the Fourier trans- 
form of the atomic density. [ll In a complementary work 
using Quantum Monte Carlo simulations, Karim Pour et 



al. interpreted the divergence of the form factors of the 
density-density and pairing correlation function as defin- 
ing a supersolid.f^l 

The purpose of this Letter is to uncover new commen- 
surability effects induced by the combination of the lat- 
tice and confinement potentials and the atomic density. 
Since the precise momentum of the ADWs in the trap 
can be controlled with the average density, it is possi- 
ble to observe effects similar to the commensurability- 
incommensurability transition appearing in crystalline 
surfaces when the density oscillations have a different 
periodicity than the underlying lattice. 0] These concepts 
have also been applied to the doped Mott transition in 
strongly correlated Fermi systems. [1, [13] We will show 
that for local densities close to half-filling there appear 
commensurate and incommensurate phases in different 
parts of the lattice. Similar local quantum criticality is- 
sues have been seen before in the repulsive case. ll| In the 
commensurate phase that arises in sectors of the lattice 
where the density is close to half-filling, the amplitude of 
the ADWs is enhanced due to local umklapp scattering 
terms. The incommensurate phase close to half- filling is 
characterized by the appearance of a new length scale, 
similar to a beating length, due to the interplay between 
the periodicity of the lattice and the periodicity of the 
ADW that results in nodes and amplitude modulation of 
the latter. 

Similar to previous works, we consider the Hamiltonian 



h.c. 



+VJ:^ {t-L/2 + Df 



(1) 



where t is the hopping, cr is a pseudospin-1/2 degree of 
freedom, U is the interaction (always attractive U < 0), 
and V is the strength of the confinement. The creation, 
destruction and number operators are the usual ones at 
each site i of the lattice. The total length L of the sys- 
tem is chosen such that the density is smooth, going to 
zero in the edges. All energies are expressed in units of t 
{t — 1). The parameter D measures the displacement of 
the center of the harmonic confinement potential with re- 
spect to the position of a lattice site. It can vary from 0.0 
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FIG. 1: DMRG and HF results for L = 180, V = 0.0005, D = 
0.0, U — —1.0, and N = 80,90. As the graph is symmetric 
only the left part is shown. The inset shows a zoom of the 
central region for A*' = 90 comparing DMRG results with 
D = 0.0 and D = 0.5. 



(the center of the harmonic well coincides with a lattice 
site) to 0.5 (the center of the well is exactly between two 
lattice sites). We will study only the unpolarized case. 

We study the ground state properties of the Hamil- 
tonian ([1]) with the DMRG algorithm[3l that provides 
very accurate numerical results. Due to the breaking 
of the translational symmetry by the trapping potential 
some modifications of the DMRG procedure are required. 
Here we follow the same procedure used before for disor- 
der potentials 1J| and umklapp scattering induced phase 
transitions. [15 In order to obtain enough accuracy for 
the largest size systems we needed to keep a maximum 
of 1400 states in each iteration. 

As a complementary tool we use the approximate 
Hartree-Fock (HF) method. It will allow us to study 
system sizes significantly beyond the limits of a DMRG 
calculation. The HF approximation describes with high 
quantitative precision the ADWs for < 1 as we have 
checked for systems up to L = 250 by comparison with 
DMRG results. For values of U between U = —1 and 
U = —2 the HF description is qualitatively good for the 
ADWs but it overestimates their amplitude. For |C/| > 2 
HF deviates from the exact results due to strong pairing 
correlations not included in the approximation. 

Without the trapping potential bosonization predicts 
an ADW with wavenumber fcADW — ^kp [5|] ■ Adding con- 
finement the wave number of the ADW is modulated by 
the density profile in the trap, but it can still be written 
as a function of an effective Fermi wavenumber fep'^J^] 



fcADW = 2fcf, fcf = ™/2, 



(2) 



an ADW with a periodicity of two lattice sites. 

Figure [1] shows DMRG results for the site density of 
a system with L = 180, U = -1, V = 0.0005, D = 0.0 
and two values of the number of atoms N = N-j + N^. 
iV = 80 is represented by open squares and = 90 by 
open circles. We also show for comparison HF results 
for the same system. The agreement between DMRG 
and HF is remarkable for t/ = — 1 . As can be seen in the 
figure, the local density is close to half-filling in the center 
of the trap. As a consequence, an ADW with a long 
range modulation of the amplitudes develops. The effect 
is more pronounced in the case with N = 90, displaying 
a nice commensurate ADW with large amplitude. The 
density in the center of the trap is below half-filling for 
A'' = 80, however the ADW can still be seen. 

The appearance of lobes and nodes in the ADW de- 
pends on small deviations of kADW from the lattice mo- 
mentum TT. These features are expected to be more pro- 
nounced for larger lattices. Therefore, we resort to the 
HF approximation to treat lattice sizes that are not at- 
tainable with present DMRG codes. In order to relate 
intermediate and large lattice systems we use the scal- 
ing NV'^ = constant while A^ ^ oo and ^ ^ 0. [i3| 
This scaling keeps constant the density in the central re- 
gion of the trap. In Fig. [5] we show HF results for the 
density in the central region for systems with L = 1000, 
V = 0.000012, U = —1 and different number of atoms, 
to illustrate the kind of ADWs appearing in larger size 
systems as a function of the number of atoms. Of par- 
ticular interest is the case of N = 560 where we can see 
the suppression of the ADW for a displacement D = 0.5 
of the trap. 

The distance between nodes 6 in the ADW is related 
to the distance it takes for the system to realize it is not 
commensurate with the underlying lattice. If the local 
average density is close to half-filling fi = 1. ± e, we have 
from ([2]) kADW — 7r(l ± e). 5 can be from the distance 
between nodes for a wave cos[(7r(l ± e)x] taken at integer 
values of the variable x, leading to 6 = e^^. 

Figure [3] shows HF and DMRG results of the distance 
between nodes 6 vs. the average local density in the cen- 
ter of the trap. The continuous line represents the for- 
mula 6 — e^^. Due to finite size effects, slight deviations 
could be seen close to half-filling. However, the overall 
agreement is excellent. 

The enhancement of the ADW amplitude close to half- 
filling can be explained in terms of the additional umk- 
lapp interaction. In the bosonization language these ef- 
fects are treated by adding a term to the umklapp part 
of the Hamiltonian. [5| 

Hs — 93 J dx cos [yS(j){x) — 5xj , (3) 



where n is the average density around the center of the 
trap. For n = 1.0 (half-filling) fcADW = and we have 



where is the coupling constant of the umklapp scat- 
tering terms, proportional to U in our model, and (t>{x) 



3 




A 0.93 

a 

V 0.9 



1 ' 1 

^ •-•N = 70 .•JjJ^fflSilulllS 
5 - flf 


1 


^ ii'ii ^ 


0<mmh/ \ \ 

I m • N = 86 ^I'flwS^^'^^'^^ 

5- si 

oLi- — — ^ — 







FIG. 2: HF local density as a function of the site index for a 
system of L = 1000 sites and different values of N and D. 
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FIG. 3: Distance between nodes (5) as a function of the av- 
erage value of the density in the center of the trap for HF 
calculations with L = 1000, V = 0.000012, and U = -1.0, 
and DMRG results with L = 100, V = 0.001, and U = -4.0. 
Inset; Amplitude of the ADW Aadw in the center of the trap 
for different values of the interaction strength with the cen- 
tral density at half-filling, L = 100, and V = 0.001. DMRG 
results (black dots), HF results for interaction strenghts up 
to U = —1.5 (diamonds). 



is the bosonized field. The additional lowering of the en- 
ergy due to the umklapp term is eflFective even for local 
densities not precisely at half-filling. This phenomenon is 
robust close to half-filling, giving rise to the enhancement 
of the amplitude and range modulation of the ADW. The 
amplitude of the ADW increases with the value of the 
attractive interaction as shown in the inset of Fig. [3] for 
DMRG and HF calculations with L = 100 and V = 0.001 
but changing the number of particles for each value of 
the interaction to keep local half-filling in the trap cen- 
ter. The increase is very abrupt for small values of |C/|, 
saturating for \U\ ~ 10. 

Another important source of commensurability comes 
from the external trapping potential. Just at half-filling 
a phase difference in the ADW can have dramatic effects 
in the observed amplitudes at the lattices sites. The cen- 



FIG. 4: DMRG results for L = 180, D = 0.0, V = 0.0005, 
and U = —4, Top panel: 86, atoms. Bottom panel: N — 70. 
The left panels show the density and the right panels show 
the local order parameters explained in the text. 



tral part of the ADW for TV = 90 is shown in the inset of 
Fig. [1] comparing the results with D = 0.0 and D ~ 0.5. 
The strong suppression of the ADW can be clearly seen. 
Even though the properties of the Luther-Emery liquid 
are not modified by the displacement of the external po- 
tential, the actual position of the lattice sites determine 
the amplitudes of the ADW. This property was not seen 
previously in the literature [l6|, because the systems stud- 
ied were too small to notice the difference but can have 
big experimental consequences for the observation of the 
ADWs. Figure [2] also shows these differences between 
D = 0.0 and D = 0.5 for a larger system. 

To proceed further with the characterization of the dif- 
ferent local phases and commensurability effects we will 
make use of two local quantities, the local variance of 
the density k,; = {nf^ — {rii)^ and the order parameter 
OcADw{i) for a local commensurate (pinned) ADW de- 
fined as: 



CADW 



1, 



(4) 



While the local variance Ki measures the fraction of spin 
paired particles in site i, the commensurate ADW order 
parameter OcADw{i) measures the local amplitude of 
the ADW between site i and its near neighbors. 

In Fig. H we show DMRG results for L = 180, ^ = 
0.0005, C/ = -4, TV = 70, and N = 86. The amphtude 
of the ADW is much higher than in the examples with 
U = — 1. The density in the case with = 70 is very 
close to half-filling in the center and the amplitude of 
the ADW very large there. The local umklapp scattering 
term effectively increases the attraction between fermions 
close to half-filling. The other case has densities larger 
than 71 = 1.0 in the central region. Nodes in the density 
due to the incommensurate ADWs appear. The ampli- 
tude of the ADWs is larger in the lobe with local density 
close to half-filling and so it is the distance between nodes 
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FIG. 5: Structure factor S{q) for different DMRG examples 
with L — 100 and V = 0.001. In the inset we show the value 
of the height of the peak at 5 = tt as a function of D for the 
example with TV = 48. 



there. In the right panels we show the value of Hi and 
OcADw{i) for both examples. The results are close to 
a plateau in the central region with a large Ki = 0.8. 
A more careful examination reveals that the maximum 
of Ki occurs in the regions close to half-filling, precisely 
where OcADwi'i) has a maximum. 

Lobes and nodes in the ADWs can be observed through 
elastic light-scattering diffraction experiments, by means 
of the Fraunhofer structure factor f3|: 



^ cxp {-iqj)nj 



(5) 



Fig.Oshows the structure factor for L — 100, V — 0.001, 
and U — —4. We show three cases with = 60 D = 0.0, 
TV = 48 D = 0.0, and Af = 48 D = 0.5. largest peak 
at g = TT corresponds to the case A^ = 48 and D = 0.0, 
exactly at half-filling. This is due to the large amplitude 
commensurate ADW. For the same number of particles 
and D — 0.5, the ADW disappears and S{q) is suppressed 
in the region around q — Above half-filling (A^ = 
60) a small peak shows up at g = tt, while most of the 
intensity concentrates at the peaks q„iax = ti""- and the 
corresponding symmetric peak. In the inset we show the 
value of ^(Tr) as a function of D to demonstrate that the 
displacement of the trap is an observable effect. 

In conclusion we have shown that commensurate and 
incommensurate phases can be seen in cold fermions 
trapped in optical lattices with attractive interaction. 
The two phases can appear spatially separated depend- 
ing on the density and confinement of the atoms in the 
lattice. The amplitude of the ADW increases with the 
strength of the interaction. The local density variance 
Ki shows that the commensurate phase is characterized 
by a maximum number of local spin-paired atoms. The 
structure factor studied in different examples reveals that 
the commensurate phase is characterized by a large peak 



exactly at fc = tt, while for phase coexistence the inten- 
sity is distributed in several peaks. The displacement 
of the trap in relation to the lattice is also a very im- 
portant parameter for the experiment. Depending on it, 
the amplitude of the ADWs and the q — tt peak can be 
suppressed. The width of the peaks is inversely propor- 
tional to the number of wells in the optical lattice. This 
parameter is the main limitation for the experimental de- 
tection of this commensurability effects. To resolve the 
two peaks appearing due to the incommensurability of 
ADWs with densities n = 1.0 ± e one needs widths of the 
peaks smaller than F < ne. According to our numeri- 
cal simulations, for harmonic confinement, the intrinsic 
width of the peak scales with size as Tint ~ 10/L. Taking 
into account that the width induced by the experimen- 
tal set-up Texp will convolute with the intrinsic width we 

shall need a system size of L > 10/^4e^/7r^ — F^^p to re- 
solve the two peaks. Currently available one-dimensional 
optical lattices of 100 sites should be enough to detect the 
results represented in Fig. [5] if the experimental width is 
less than 0.24 in units of the density[3|]. 
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